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1 We briefly review the Thomas-Fermi statistical model of atoms in the classical 

H ' non-relativistic formulation and in the generalised finite-nucleus relativistic formu- 

^ ' lation. We then discuss the classical generalisation of the model to finite tempera- 

ff^ ' tures in the non-relativistic approximation and present a new relativistic model at 

^ finite temperatures, investigating how to recover the existing theory in the limit 

, of low temperatures. This work is intended to be a propedeutical study for the 

• evaluation of equilibrium configurations of relativistic "hot" white dwarfs. 

o : 

' 1 Introduction 

The Thomas-Fermi statisticaljjnodel of the atom has been extensively used, since 
I ' its early formulation in 1927 Id'El, to evaluate the equation of state of compressed 

2 , matter and in general in the description of the electronic structure^ of atoms and 

' solids under different physical conditicjiia El. In our previous work □ we applied a 

, relativistic generalisation of this model Era to the evaluation of the equation of state 

^ ' and the equilibrium configurations of cold white dwarfs, whereas "cold" means that 

. ^ , the temperature of the star is much smaller than the Fermi energy of the degenerate 

' electron gas. We are now interested in the finite temperature effects on the equation 

of state and we thus derive a temperature dependent formulation of the relativistic 
Thomas-Fermi model. The paper is organised as follows: first (section 2) we briefly 
describe the classical Thomas-Fermi model, i.e. the non-relativistic and completely 
degenerate case. Then (section 3) we'll review the finite temperature formulation 
of this problem , giving the corresponding exact equation and its approximated 
version. We'll thus pass to the relativistic theory in the completely degenerate case 
(section 4) and finally to the new relativistic model at finite temperature (section 
5). 



2 Classical Thomas-Fermi model 

We first consider the simple Thomas-Fermi model. Let us consider the spherically 
symmetric problem of a nucleus with Z protons and A nucleons interacting with a 
fully degenerate gas of Z electrons. 



procmgO: submitted to World Scientific on February 1, 2008 



The fundamental equation of electrostatics for this problem is 

AV{r) = 47rene(r) (1) 

where V{r) is the electrostatic potential and ne{r) is the number density of elec- 
trons. 

The electrostatic potentialis related to the Fermi momentum (and thus to the 
number density of electrons) by the equilibrium condition (see a) 

p%{r) /2m - eV{r) = const = Ef (2) 

where the name Ep stands for Fermi-Thomas chemical potential or Fermi Energy 
of the electrons. 

To put the equation |l| in adimensional form we introduce the new function <i>(r), 
related to the coulomb potential by 

<^{r)^V{r)+EF/e (3) 

and the corresponding adimensional function x, implicitly defined by 

$(r) = ^ (4) 
r 

Furthermore we introduce the new independent variable x, related to the radius 
r by the relation r — bx, where 

It is easy to show that in this case one can write eq.^ in the form 



(6) 



which is the classical adimensional form of the Thomas- Fermi equation. 

The first initial condition for this equation follows from the request that ap- 
proaching the nucleus one gets the ordinary Coulomb potential 

X(0) = 1 (7) 
The second condition comes from the normalisation condition 

N = Anrier'^dr (8) 



Jo 

which gives 

N = Z[xox'{xo)-xM + l] (9) 
and for neutral atoms {N = Z) 

xox'i^o) = xi^o) (10) 
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3 Temperature Dependent Non-relativistic Model 



We introduce now the finite temperature effects in the modei. This was aiready 
done in 1940 by Marshack and Bethe Q through a perturbation treatment, while the 
full adimeiisional equation is discussed in a successive work of Feynman, Metropolis 
and Teller u. 

The density of an electron gas at temperature T can be written as 

ne = - 1-/ — de^- -:^{KTf'^h[4—] 11 

TT-^K" Jo e^+1 ^2;j3 V ) ly^j,) V ) 

where 

Ux) = / — ^^^dy (12) 

Using the same adimensional variables introduced in the previous section and 
introducing the temperature parameter r as 

r=^KT (13) 
we can rewrite the electron density in adimensional form 



h (—) (14) 



We can thus express the electrostatic equation in the following adimensional 
form 

This equation is formally different from the one obtained in the case of complete 
degeneracy, nevertheless it can be easily shown that if one develops the integral 
which appears in eq. ^ for small temperatures (see appendix for details) one gets 
the following formula at the first order 



1 



8 



(16) 



dx^ x^/"^ 
where we neglect terms of the order 0{t'^). 

4 Relativistic model at T=0 

When considering a relativistic extension of the model the finite size of the nucleus 
must be taken in account to avoid the central singularity (see e.g. □). In this case we 
have thus to introduce a new term into the fundamental equation of electrostatics, 
representing the positive distribution of charges in the interior of the nucleus np{r) 



/S.V{r) = 47rene(r) - 4:TTenp{r) (17) 
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The quantities V(r) and n(r) are in this case related by 



c\J p^p + m?c^ — eV{r) — const = Ep 
Using eq.|| it is possible to put eq. |lj in the form 



-<i>^ + 2me$ 



which, using y,becomes 

Pf = "imc 



Z 

Zcr 



where 



2/3 



1/2 



1/2 



z 

Zf-r 



4/3 



1/2 



Zrr=\—\ ( ^ 



3/2 



2462.4 



(18) 



(19) 



(20) 



(21) 



Remembering the relation between the Fermi momentum and the number den- 
sity of a fermion gas 



Pi 



we obtain the following expression 



(- 



3/2 



z 



4/3 



3/2 



(22) 



(23) 



We can also express the second term of the right-hand side of eq.|i| in terms of 
adimensional quantities: we assume here an homogeneous spherical nucleus, with 
a radius given by the approximate formula 



rnuc = 1.2v4i/3 fm 
The number density of protons is therefore 

3Z 



Finally we can write eq.O in the form 



(fx ^ 



Z 

Zcr 



-1 3/2 



X) 



3a; 



■e(x. 



X) 



(24) 



(25) 



(26) 



where Xnnc is the adimensional size of the nucleus (r„„c = bxnuc)- 

Equation |6| is what we call ^Generalised adimensional Fermi-Thomas equation\ 
The first initial condition for this equation follows from the fact that x 

and therefore x ~^ 0, and so 



x(o) = 



(27) 
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The second condition comes from the normaUsation condition 

3/2 

X dx (28) 



^0 



X 



1/2 



1 



^ \ 4/3 



with To — bxQ atom size. Developing this formula we have 



N = Z I xx" dx + ^— / x^ dx + Z I xx" dx (29) 



Xo 



3Z 

3 

which gives again the relation 

7V-Z[xox'(2:o)-x(xo) + l] (30) 

Note that the physical quantities such as the coulomb potential and the density 
of electrons do not show any singularity in the center, neither on the border of the 
nucleus, being dependent just on the function x his first derivative. The only 
discontinuity appears in the second derivative of x due to our rough assumption of 
homogeneous spherical nucleus. 

It is also evident that the scaling properties of the classical Thomas-Fermi equa- 
tion are lost in this case, where one has to integrate the adimensional equation 
separately for each value of Z and different values of xq, i.e. for different states of 
compression. 

In figg.l and 2 we show an example of the applications of this model to the 
||tudy of equilibrium configurations of cold White Dwarfs. For more details see ref. 

5 Temperature Dependent Relativistic Model 

We consider now the complete problem of a relativistic and degenerate gas of elec- 
trons at temperature T surrounding a positively charged nucleus. 
The number density of such a gas is 



TT^ichf J^,. + 1 ' " 7r2(cn)3 2 \kt) 



(31) 



where 



and a is the fixed parameter a — mc^ / (KT) . 

Introducing adimensional variables and the temperature parameter r we can 
write now 



1 Z^ 



2\ 3 



(— ) (33) 



'■2 

TX 



and inserting the previously defined quantity Z^r we find 



47rfe3 \Z. 



TX 
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Figure 1. Equilibrium configurations curve for Oxygen WD in the M/Mq — pc plane, obtained 
in Newtonian theory and General relativity. For comparispji we show the results obtained with 
Chandrasekhar model .Reproduced from Bertone & RuffiniQ. 



The final and more general expression of the Thomas-Fermi equation is thus 
found to be 




This formula is the main result of this paper and will be the point of depar- 
ture for the evaluation of the equation of state of compressed matter in extreme 
conditions of high temperatures in a relativistic regime. 

We will present elsewhere the numerical integration of this equation and 
the corresponding equilibrium configurations of hot relativistic white dwarfs. 

6 Conclusions 

We have first discussed the classical Thomas-Fermi method and the corresponding 
generalisation to finite temperature or to the relativistic regime. A finite nucleus 
treatment including both effects has then been presented. This work must be 
considered as propedeutical to the evaluation of the equilibrium configurations of 
relativistic hot white dwarfs. The numerical integration of the highly non-linear 
equation obtained and the application to the equation of state of white dwarfs 
matter will be presented elsewhere. 
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Figure 2. Equilibrium configurations curve for Iron WD in the M/Mq — pc plane, obtained in 
Newtonian theory (the General relativistic curve is pratically superposed at the newtonian one at 
these low densities). For comparison me show the results obtained by Chandrasekhar and Salpeter 
.Reproduced from Bertone & RufBnO. 



Appendix 

Following Landau and Lifshitz 1^ we recall here how to approximate for low tem- 
peratures integrals of the form 

■^^^^ -de (36) 



e KT + I 

appearing in the statistical treatment of fermions. 

Without going into details, we just recall that using the fact that fi/T >> 1 
one can put the former integral in the form 

7o + 1 3 7o + 1 

(37) 

These integrals can be evaluated as 

/ ^-T^^=(l~2i--)r(x)^- (38) 

^ n— 1 
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We thus obtain our approximate result 

/ = £ f{e)de + ^{KTffifj,) + I^(i^T) + ••. (39) 
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